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This paper studies the design of a decentralized multiuser multi-antenna (MIMO) system for spectrum sharing 
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' over a fixed narrow band, where the coexisting users independently update their transmit covariance matrices for 

individual transmit-rate maximization via an iterative manner. This design problem was usually investigated in 
the literature by assuming that each user treats the co-channel interference from all the other users as additional 
(colored) noise at the receiver, i.e., the conventional single-user decoder (SUD) is applied. This paper proposes a 
new decoding method for the decentralized multiuser MIMO system, whereby each user opportunistically cancels 
the co-channel interference from some or all of the other users via applying multiuser detection techniques, thus 
termed opportunistic multiuser detection (OMD). This paper studies the optimal transmit covariance design for 
^ ■ users' iterative maximization of individual transmit rates with the proposed OMD, and demonstrates the resulting 

capacity gains in decentralized multiuser MIMO systems against the conventional SUD. 
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O ' I. Introduction 

O 

The Gaussian interference channel is a basic mathematical model that characterizes many real-life 



^ communication systems with multiple uncoordinated users sharing a common spectrum to transmit 



independent information at the same time, such as the digital subscriber line (DSL) network [1], the 
ad-hoc wireless network [2], and the newly emerging cognitive radio (CR) wireless network [3]. From 
an information-theoretical perspective, the capacity region of the Gaussian interference channel, which 
constitutes all the simultaneously achievable rates of the users in the system, is still unknown in general 
[4], while significant progresses have recently been made on approaching this limit [5], [6]. Capacity- 
approaching techniques usually require certain cooperations among distributed users for their encoding 
and decoding. A more pragmatic approach that leads to suboptimal achievable rates of the users in the 
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Gaussian interference channel is to restrict the system to operate in a decentralized manner [7], i.e., 
allowing only single-user encoding and decoding by treating the co-channel interference from the other 
users as additional Gaussian noise at each user's receiver. In such a context, decentralized algorithms 
for users to allocate their transmit resources such as the power, bit-rate, bandwidth, and antenna beam 
to optimize individual transmission performance and yet to ensure certain fairness among all the users, 
become most important. 

This paper focuses on a multiuser multiple-input multiple-output (MU-MIMO) wireless system, where 
multiple distributed links, each equipped with multiple transmit and/or receive antennas, share a common 
narrow band for transmission in a fully decentralized manner. In such a scenario, the system design reduces 
to finding a set of transmit co variance matrices for the users subject to their co-channel interference result- 
ing from their simultaneous and uncoordinated transmissions. This design problem has been investigated 
in a vast number of prior works in the literature, e.g., [8]-[16], by treating the co-channel interference as 
additional colored noise at each user's receiver, i.e., the conventional single-user decoder (SUD) for the 
classic point-to-point MIMO channel is applied. In [8], the authors proposed an algorithm, which is in 
spirit analogous to the iterative water-filling (IWF) algorithm in [7], for each distributed MIMO link to 
iteratively update transmit covariance matrix to maximize individual transmit rate. Distributed iterative 
beamforming (the rank of transmit covariance matrix is restricted to be one) algorithms were also studied 
in [9] for transmit sum-power minimization given individual user's quality of service (QoS) constraint 
in terms of the received signal-to-interference-plus-noise ratio (SINR). The throughput of decentralized 
MU-MIMO systems has been further analyzed in [10] and [11] for the cases of fading channels and large- 
size systems, respectively. In [12], [13], centralized strategies were proposed where all users' transmit 
covariance matrices are jointly searched to maximize their sum-rate, and numerical algorithms were also 
proposed to converge to a local sum-rate maxima. Analyzing the decentralized MU-MIMO system via a 
game theoretical approach has recently been done in [14]- [16]. 

The cited papers on decentralized/centralized designs for the Gaussian MIMO interference channel have 
all adopted the SUD at each user's receiver, whereas during the past decade multiuser detection techniques 
(see, e.g., [17] and references therein) have been thoroughly investigated in the literature, and have been 
proven in realistic multiuser/MIMO systems to be able to provide substantial performance gains over 
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the conventional SUD. This motivates our work's investigation of the following question: Considering 
a decentralized MU-MIMO system where the users iteratively adapt their transmit covariance matrices 
for individual rate maximization, "Is applying multiuser detection at each user's receiver able to enhance 
the system throughput over the conventional SUD?" Note that because of the randomness of channels 
among the users, as well as their independent rate assignments, at one particular user's receiver, multiuser 
detection can be used to cancel the co-channel interference from some/all of its coexisting users only 
when their received signals are jointly decodable with this particular user's own received signal. Thus, 
we refer to this decoding method as opportunistic multiuser detection (OMD). Also note that the OMD 
in the context of the decentralized MU-MIMO system is analogous to the "successive group decoder 
(SGD)" in the fading multiple-access channel (MAC) with unknown channel state information (CSI) at 
the user transmitters (see, e.g., [18] and references therein). With the proposed OMD, this paper derives 
the optimal transmit covariance matrix for user's individual transmit-rate maximization at each iteration 
of transmit adaptation. By simulation, this paper demonstrates the throughput gains of the converged 
users' transmit covariance matrices with the proposed OMD over the conventional SUD. 

The rest of this paper is organized as follows. Section |Il] presents the system model of the decentralized 
MU-MIMO system. Section [nl] studies the optimal design of user transmit covariance matrix with the 
proposed OMD for the special case with two users in the system. Section |W] generalizes the results to 
the case of more than two users. Section |V] provides the simulation results to demonstrate the throughput 
gains with the proposed OMD over the SUD. Finally, Section |VI] concludes the paper. 

Notation: Scalars are denoted by lower-case letters, e.g., x, and bold-face lower-case letters are used 
for vectors, e.g., x, and bold-face upper-case letters for matrices, e.g., X. In addition, tr(5'), \S\, S~^, 
and 52 denote the trace, determinant, inverse, and square-root of a square matrix S, respectively, and 
S ^ means that S is a positive semi-definite matrix [19]. For an arbitrary- size matrix M, denotes 
the conjugate transpose of M. diag(xi, . . . , xm) denotes a M x M diagonal matrix with xi, . . . , xm as 
its diagonal elements. / and denote the identity matrix and the all-zero vector, respectively. E[-] denotes 
the statistical expectation. The distribution of a circular symmetric complex Gaussian (CSCG) random 
vector with mean x and covariance matrix S is denoted by CJ\f{x, S), and ~ stands for "distributed as". 
C^^y denotes the space of xxy matrices with complex-valued elements. ma,x(x, y) and min(a;, y) denote 
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the maximum and minimum between two real numbers, x and y, respectively, and = max(x, 0). 

II. System Model 

This paper considers a distributed MU-MIMO system where K users transmit independent information 
to their corresponding receivers simultaneously over a common narrow band. Each user is equipped with 
multiple transmit and/or receiver antennas, while for user k, k — 1, . . . , K, Nf. and Mj. denote the number 
of its transmit and receive antennas, respectively. For the time being, it is assumed that perfect time and 
frequency synchronization with reference to a common clock system have been established for all the 
users in the system prior to data transmission. We also assume a block-fading model for all the channels 
involved in the system, and a block-based transmission for all the users over each particular channel 
fading state. Since the proposed study applies to any channel fading state, for brevity we drop the index 
of fading state here. The discrete-time baseband signal for the A;th user transmission is given by 

K 

Vk = HkkXk + ^ HjkXj + Zk (1) 

j=l,j^k 

where Xk e C^*'^^ and y,, e c^'^^^ are the transmitted and received signal vectors for user k, respectively, 
A; e {1,. ..,K}; Hkk G C*^''-^^'= denotes the direct-link channel matrix for user k, while Hjk G C^*"^^^ 
denotes the cross-link channel matrix from user j to user fc, j e {1, . . . , K}, j ^ k; and Zk G C^*^^ is 
the received noise vector of user k. 

Without loss of generality, it is assumed that Zk ~ CJ^{0, 1),\/k e {1,...,K}, and all Zk's are 
independent. We consider a decentralized multiuser system where the K users independently encode 
their transmitted messages and thus x^'s are independent over k. Since this paper is interested in 
the information-theoretic limit of each Gaussian MIMO channel involved, it is assumed that Xk ~ 
CM{0, Sk),yk e {1, . . . , K}, where Sk — 'E[xkxj^] is the transmit covariance matrix for user k. 

This paper considers a similar decentralized operation protocol as in [7], [8], [14]-[16], whereby the 
users in the system take turns to update their transmit covariance matrices for individual rate maximization, 
with all the other users' transmit covariance matrices being fixed, until all users' transmit covariance 
matrices and their transmit rates get converged. We consider two types of decoding methods at each 
user's receiver. One is the conventional SUD, which has been applied in the above cited papers, where 
the kth user decodes its desired message by treating the co-channel interference from all the other users. 
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j 7^ k, as additional colored Gaussian noise ~ CAf{0, X] 



.K 



HjkSjHjf^). The other decoding method 



is the newly proposed OMD, whereby each user opportunistically applies multiuser detection to decode 
some/all of its coexisting users' messages so as to cancel their resulted interference, provided that these 
messages are jointly decodable with this user's own message. In practice, each user in the system is 
usually interfered with by all the other users, while due to location-dependent shadowing/fading, only a 
small group of coexisting users who are closest to one particular user and thus correspond to the strongest 
cross-link channels to this user, will contribute the most to this user's received co-channel interference. 
As a result, this user can effectively estimate the transmit rates as well as the cross-link channels of these 
"strong" interference users, and employ the proposed OMD to suppress their interference at the receiver. 
Note that the use of OMD instead of SUD still maintains the fully decentralized property of the existing 
IWF-like operation protocols given in [7], [8], [14]-[16]. 



In this section, we present the problem formulation as well as the solution to determine the optimal 
transmit covariance matrix of each user for individual transmit-rate maximization, when the proposed 
OMD is employed. For the purpose of exposition, we consider the special case where only two users 
exist in the system. We will address the general case with more than two users in Section |IVl For brevity, 
only user I's transmit adaptation is addressed here, while the developed results apply similarly to user 2. 

A. Problem Formulation 

Note that at one particular iteration of user 1 to update its transmission, user 2's transmit covariance 
matrix, S2, and transmit rate, denoted by r2, are both fixed values. For a given transmit covariance matrix 
of user 1, Si, the resultant maximum transmit rate of user 1 can be expressed as 



III. Transmit Covariance Optimization: The Two-User Case 




(2) 



where 



R2 



.{a) 



log |I + (/ + HiiSiH^i)-'H2iS2H. 



H 
21 



(3) 



i?2 



\og\l + H21S2H. 



(4) 



The above result is illustrated in the following three cases corresponding to the three expressions of 
ri in Q from top to bottom. 

• Strong Interference Case: In this case, the received signal from user 2 is decodable at user I's 
receiver with the conventional SUD, by treating user I's signal as colored Gaussian noise. This 
is feasible since r2 < b!^' given in ([3]). After decoding user'2 message and thereby canceling its 
associated interference, user 1 can decode its own message with a maximum rate equal to its own 
channel capacity. The above decoding method is known as successive decoding (SD) for the standard 
Gaussian MAC [20]. 

• Moderate Interference Case: In this case, r2 > R^2^ and thus the received signal from user 2 is not 
directly decodable by the SUD. However, since r2 < i?2^^ given in (U), it is still feasible for user 1 
to apply joint decoding (JD) [20] to decode both users' messages!!] In this case, the rate pair of the 
two users should lie on the 45-degree segment of the corresponding MAC capacity region boundary 
[20], i.e., ri + r2 = log |/ + HnSiH^, + H^iS^H^^^ 

. Weak Interference Case: In this case, r2 > i?2^\ i.e., the received signal from user 2 is not decodable 
even without the presence of user I's signal. As such, user I's receiver has the only option of treating 
user 2's signal as colored Gaussian noise and applying the conventional SUD to directly decode user 
I's message, the same as that in the existing IWF-like algorithms (see, e.g., [8], [14]-[16]). 
In the above decoding method, multiuser detection is applied in both cases of strong and moderate 
interferences when r2 < \ but not in the case of weak interference when r2 > i?2^'*. Thus, user I's 
receiver opportunistically applies multiuser detection to decode user 2's message, either successively (SD) 
or jointly (JD) with its own message. We thus refer to this decoding method as opportunistic multiuser 
detection (OMD). From ([3]) and (H)), it follows that i?2°'' < ^2^^- Further more, it is easy to verify that ri 
given in Q with the proposed OMD is in general larger than the achievable rate with the conventional 
SUD (given by the third expression of ri in Q independent of r2), for any given set of Si, ^2, and r2. 
With ri(S'i) given in (O for a fixed Si, we can further maximize user I's transmit rate by searching 

'Note that SD can also be applied in this case to achieve the same rate for user 1 as JD, if SD is deployed jointly with the "time sharing" 
[20] or "rate splitting" [21] encoding technique at user I's transmitter Since these techniques require certain cooperations between users, 
they might not be suitable for the fully decentralized multiuser system considered in this paper. 



over Si. Let Pi denote the transmit power constraint of user 1. This problem can be expressed as 

(PI) max ri(5i) 
Si 

s.t. tr(5i) < ^ 

where ri(S'i) is given in The optimal solution of Si in (PI) and the corresponding maximum transmit 
rate of user 1 are denoted by 5^^° and rP^°, respectively. 

B. Proposed Solution 

In this subsection, we study the solution of (PI) for the optimal transmit covariance matrix of user 1, 
when the proposed OMD is deployed at user 1' receiver. Note that although the constraints of (PI) are 
convex, its objective function is not necessarily concave due to the fact that R'^^ given in ([3]) is neither 
convex nor concave function of Si. As a result, (PI) seems to be non-convex at a first glance. In fact, 
(PI) is a convex optimization problem after being transformed into a convex form, as will be shown 
in this subsection. In the following, we will study the solution of (PI) for two cases: r2 > R'^^ and 
r2 < Rf\ for which the SUD and the multiuser decoding (MD) (in the form of either SD or JD) should 
be used to achieve ri(S'i) given in (O, respectively. 

1) r2> R'il- In this case, the SUD should be applied. Note that R^^ is a constant unrelated to Si. 
Thus, the optimal Si that maximizes the third expression of ri{Si) in Q has the following structure 
[20]: 

^SUD _ YJ^yH 

where V E ^^ixTi ^^^^i Ti = min(A''i,Mi) is obtained from the singular-value decomposition (SVD) 
of the equivalent channel of user 1 (after the noise whitening) expressed as 

(/ + H2iS2Hg)-"^Hii = C/SV^ (6) 

with U G C^^i""^!, S = diag(cri, . . . , (TtJ, (Ji > 0, i = 1, . . . , Ti, and A = diag(pi, . . . with pi's 
obtained from the standard water-filling solution [20]: 

Pi=(fi-\] , i = l,...,Ti, (7) 



with fi being a constant to make YliLiPi = Pi- The maximum rate of user 1 then becomes 

^'^°-Elog(l + ^.^P.)- (8) 



1=1 



2) r2< Rf: In this case, the MD in the form of either SD or JD should be used. In order to overcome 
the non-concavity of ri(S'i) given in ^ due to R^-f', we re-express the first two expressions of ri(S'i) 
in © as 

rf °(5i) = min (log |/ + //ii5iiffi| , log |/ + /fii5i/ffi + H2iS2Hg\ - rs) . (9) 
Thus, the maximum achievable rate of user 1 can be obtained as 

rf°= max rf°(5i). (10) 

S'i:tr{S'i)<Pi,>S'ihO 

The maximization problem in (flOl ) can be explicitly written as 
(P2) max ri 

ri, Si 

s.t. n<\og\l + HnSiH^,\ (11) 

n <log\l + HnSiH"^^ + H^iS^Hg] ~ r2 (12) 

ri > 0,tr(5i) < Pi,5i ^ 0. (13) 

The optimal solution of ri in (P2) will be r^^. Note that (P2) is a convex optimization problem since 
its constraints specify a convex set of {ri,Si). To solve (P2), we apply the standard Lagrange duality 
method [19]. First, we introduce two non-negative dual variables, /ii and /i2, associated with the two rate 
constraints (fTTl) and (fT2l) . respectively, and write the associated Lagrangian of (P2) as 

£(ri,5i,/ii,^2) =n -/ii (ri - \og\l + HuSiHg\) 

- fi2 (ri - log |/ + HnSiHg + /fai^s/ffi] + rs) (14) 
By reordering the terms in (fT4l) . we obtain 

/:(ri, 5i,/ii,/i2) =(1 - yUi - /i2)ri + ^1 log |/ + HnSiHf^\ 

+ ^2log |/ + HnSiHg + H2iS2Hg\ + /iars- (15) 
The Lagrange dual function of (P2) is then defined as 

g{jj,i,fi2)= max C{ri, Si, jJi, ^2) (16) 
{ri,Si)eA 



where the set A specifies the remaining constraints of (P2) given in (fT3l) . The dual problem of (P2), of 
which the optimal value is the same as that of (P2)o is defined as 

(P2-D) rnin a(/ii,/i2). (17) 

Mi>0,M2>0 

Let rl and denote the optimal solutions of (P2). Let /i* and ^2 denote the optimal dual solutions 
of the dual problem (P2-D). Next, we will present a key relationship between /i* and /ig as follows. 

Lemma 3.1: In problem (P2-D), the optimal solutions satisfy that + jj^ = 1- 

Proof: See Appendix HI ■ 

Given Lemma I3.1[ without loss of generality, we can replace by 1 — /ii in (fT5l) . Thus, the 
maximization problem in (fT6l) can be equivalently rewritten as (by discarding the constant term fi2r2) 

(P3) max ^ilog|/ + JfnSiiffil + (l-/ii)log|/ + Jfii5iiffi + i/2i52Jf^i| 
Si 

s.t. tr(5i) < ^ 0. (18) 

Further more, the dual problem ([17] ) now only needs to minimize (^(/ii) (since /i2 = 1 — fJ-i) over 
< /Ui < 1. Then, there are the following three cases in which f^l takes different values. 

• fil = 0: In this case, fi^ = 1. From the Karush-Kuhn-Tucker (KKT) optimality conditions [19] of 
(P2), it is known that the constraint ([TT]) is inactive while the constraint ([T2[) is active. This suggests 
that JD instead of SD is optimal. Furthermore, from (P3), with jji = fi^ = 0, it follows that S*^, 
denoted by S'l^, maximizes the sum-rate, log |/ + HnSiH^^ + H2iS2H2i\, from which we can 
show that 

5f = 5f ° (19) 

where Sf^^ is given in ©, i.e., the optimal transmit covariance matrix is the same for both cases 
of SUD and JD. However, the optimal r* in this case with JD, denoted by rf^, is equal to 

r f = rSUD + r(^) _ r2 (20) 

where rf^^ is given in ([8]). Finally, we need to check the condition under which this case holds. 
Since the constraint ([TT]) should be inactive, it follows that 

<log|/ + //n'S'fJ/fi|. (21) 

^It can be easily checked that the Slater's condition holds for (P2) and thus the duality gap for (P2) is zero [19]. 
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From (|20l) and (|2T1) . it can be shown that the case of interest holds when 

r2 > log |/ + (/ + /fn5f if fi)-i/f 21 -S-s/ffi I = R^^^ (22) 

Note that ^2"'* can also be obtained from i^g"^-* given in (|3]) by letting Si = 5^°. 
;U* = 1: In this case, = 0. From the KKT optimality conditions of (P2), it is known that the 
constraint (fTTI) is active while the constraint (fT2l) is inactive. This suggests that SD instead of JD is 
optimal. Furthermore, from (P3), with ni = = 1, it follows that S\, denoted by , maximizes 
user I's own channel capacity (without the presence of user 2), log | / + HiiSiH^i\, from which 
we can easily show that [20] 

Sf = ViAiVf (23) 

where Vi E C^^^^^^ is obtained from the SVD of the direct-link channel of user 1 expressed 
as Hn = C/irVf, with G C*^^x^S Ti = diag(7i, . . . , 7tJ, % > 0, t = l,...,Ti, and 
Ai = diag(q'i, . . . , grj with g^'s obtained from the standard water-filling solution [20]: 

qi= ^ly- -^^ , i = l,...,Ti, (24) 

with u being a constant to make J2j=i Qi = Pi- The optimal in this case with SD, denoted by 
rf^, then becomes 

Ti 

rf° = J]log(l + 7fg.). (25) 

i=l 

Similarly like the previous case, we can show that this case holds when 

r2 < log |/ + (/ + iffi)-^/f2i52/ffi| = Ri"\ (26) 

At last, we have the following lemma. 

Lemma 3.2: For R^^^ defined in and R^^^ defined in it holds that i?^"^ > R^^^K 

Proof: See Appendix HIl ■ 
< /i* < 1: In this case, < /ig < 1, and from the KKT optimality conditions of (P2), it is known 
that both the constraints ([TT]) and (fT2l) are active. This suggests that rl = log [/ + HnSlHf^^, i.e., 
SD is optimal. However, the optimal solution of (P2), or that of (P3) with fii = fil, denoted 

~ SD 

by Si , in general does not have any closed-form expression, and thus needs to be obtained by 
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a numerical search. Since (P3) is convex, the interior-point method [19] can be used to efficiently 
obtain its solution for a given fii. Let S\{fii) denote the optimal solution of (P3) for a given /ii. 
Then, /i* can be efficiently found by a simple bisection search based upon the sub-gradient [19] of 
g{fii), which can be shown from ([T?] ) (with fj,2 = I — fJ-i) to be 

I + {I + HnSl{fii)H^,y' H2iS2H^i 



log 



r2- 



(27) 



~ SD 



Once fii converges to /i^, the corresponding SKfii) becomes the optimal . The optimal in 
this case with SD, denoted by fp, is then expressed as 



log 



I + HiiSi i3"^ 



(28) 



Similarly like the previous two cases and using Lemma I3.2[ we can show that this case holds when 



<r2<M"'. 



(29) 



3) Combing r2 > R2^ and r2 < Rf\' To summarize, the following theorem is obtained for the optimal 
solution of (PI). 

Theorem 3.1: For a given set of ^2 and r2 of user 2, the optimal transmit covariance matrix of user 
1 and the maximum transmit rate of user 1 with the proposed OMD are given as follows: 



OMD 



Sf, < r2 < R^2^ 



^SD 

"^1 5 



Si 



JD 



Ri"^ <r2< Rf 
R^^ <r2< ^ 



(30) 



5? 



SUD 



r2 > R 



(b) 
■2 1 



,OMD 
1 



».SD 
' 1 ' 

^SD 
' 1 ) 

».JD 
'1 ) 

^SUD 
' 1 > 



< r2 < R^2^ 
R^^ <r2< R^2^ 
R^^^ <r2< Rf 
r2>Rf. 



(31) 



The corresponding optimal decoding methods at user I's receiver are (from top to bottom) SD, SD, JD, 
and SUD, respectively. 

In Fig. [U we show in (|3T]) as a function of r2 for some fixed 82- The rate gain of rf^^ for OMD 



over rf^^ for SUD is clearly shown when r2 < R2' ■ There are three pentagon- shape capacity regions 



shown in the figure, which are Cmac('S'^°, 'S'2), Cmac{S\" , S2), and CyiAc{S'{" ■, S2), respectively, where 



SD 



~SD 
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Cmac('S'i, S2) denotes the capacity region of a two-user Gaussian MIMO-MAC with user I's and user 
2's transmitters transmitting to user I's receiver, and Si, S2 denoting the transmit covariance matrices 
of user 1 and user 2, respectively. More specifically, Cmac('S'i, 'S'2) can be expressed as [20] 



Cmac{Si, S2) = l{n, r2) : ^ Ti < log 



yj<Z{l,2}\. (32) 



Note that in Fig. [H the sold line consisting of different rate pairs of [r'^^^ , r2) constitute the boundary 



rate pairs of the aforementioned capacity regions. Also note that there is a curved part of this rate-pair 
line in the case of i^g''^ < r2 < R'^\ where rf^^ is equal to rf^ and is achievable by , which is the 
solution of problem (P3) for some given yUi, < /xi < 1. 

IV. Extension to More Than Two Users 

In this section, we extend the results obtained for the two-user MIMO system to the general MU-MIMO 
system with more than two users, i.e., K > 2. Due to the symmetry, we consider only user I's transmit 
optimization over Si to maximize transmit rate ri, with all the other users' transmit rates, r2, . . . , r^^, 
and transmit covariance matrices, ^2, • • • , Sk, being fixed. 

To apply OMD at user I's receiver, we need to first identify the group of users whose signals are 
(jointly or successively) decodable at user I's receiver without the presence of user I's own received 
signal. We thus have the following definitions: 

Definition 4.1: A set Wi, Wi C {2, . . . , K}, is called a decodable user set for user 1, if the received 
signals at user I's receiver due to the users in Ui are decodable without the presence of user I's own 
received signal, by treating the received signals from the other users in Ui as colored Gaussian noise, 
where Ui denotes the complementary set of Ui, i.e., Uif]Ui = and Ui\jUi = {2, . . . , K}. More 
specifically, the transmit rates of users in Ui must satisfy [20] 

1+ ll+J2^^kiSkH^i\ ^/f,i5,f/,^ ,\/JCUi. (33) 



< log 



Definition 4.2: A set C {2, . . . , K} is called an optimal decodable user set for user \, if is a 
decodable user set for user 1, and among all possible decodable user sets for user 1, has the largest 
size. 

Next, we have the following important proposition: 
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Proposition 4.1: The set U* is unique. Furthermore, for any decodable user set for user 1, Wi, it holds 
that^Yi CU*. 

Proof: See Appendix Hill ■ 
For conciseness, we show the algorithm to find the unique set for user 1, U*, in Appendix HVl 
From Proposition 14. 1[ it follows that the optimal decoding strategy for user I's receiver is applying 
OMD to the users in the set (it may be possible that = 0), while taking the users in the set as 
additional colored Gaussian noise. For an arbitrary set V, let |V| denote the size of V. Note that to make 
the OMD feasible, the rate of user 1, ri, and the rates of users in must be jointly in the capacity region 
of the corresponding {\V(l \ + l)-user Gaussian MIMO-MAC for a given set of user transmit covariance 
matrices and the receiver noise covariance matrix, $ = J + J2k&n* -^kiSkHj!^, which, similar to (|32l) . 
can be defined as 



Cmac(Wi) =<^ (ri, {ri}ieu*) : ^ r ^ < log 



ieJ 



,VJC{1}|JW*|. (34) 



Note that in (l34l) . the rate inequalities involving subsets jT's containing users solely from all hold due 
to the definition of Wj*. Therefore, in order to find the optimal Si for user 1 to maximize ri, with fixed 
Tj's and Si's, i = 2, . . . , K, it is sufficient to consider the following optimization problem: 



(P4) max ri 

Si,ri 



S.t. ri + ^ Ti < log 

ri > 0,tr(5i) < Pi, Si h 



1+^ ' lHiiSiHf, + J2HiiS,HfA 

V ieJ J 



yjCU* (35) 
(36) 



Problem (P4) is convex in terms of ri and Si since its constraints specify a convex set of (ri,S'i). 



Similarly like for problem (P2), we introduce a set of non-negative dual variables, /i„'s, n 



2l"i*l 
, . . . , ^ , 

each associated with one corresponding constraint in (1351) for a particular subsect J' (including JT" = 0) 
denoted by J'n, and obtain an equivalent problem for the optimization over Si for a given set of fixed 
/x„'s, which is expressed as 



(P5) max ^ Hn log 
(-J 1 



n=l 



/ + I HiiSiH^, + ^ H,iS,Hf, 



s.t. tr(5i) < Pi,5i ^ 0. 



(37) 
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It can be shown that problem (P5) is convex, and thus it can be solved via standard convex optimization 
techniques, e.g., the interior point method [19], while in general, no closed-form solution for (P5) is 
available, similar to the previous two-user case in Section Unl Let the optimal solution of (P5) be denoted 
by Sl{{n„}). Then, /i„'s can be updated towards the optimal dual solutions of (P4) via the well-known 
ellipsoid method [19] subject to an additional constraint, J2nf^n = 1 (similar to Lemma [3711 in the two- 
user case). Let the optimal solutions of ^„'s be denoted by yU*'s. The optimal solution of Si for (P4) 
with OMD is then obtained as 

^OMD ^ si{{n*J), and the corresponding maximum achievable rate of 
user 1, rP^°, can be obtained from any active constraint in (l35l) with equality. The optimal decoding 
orders/decoding methods for the users in prior to decoding user I's message can be obtained according 
to the optimal non-zero dual solutions, //* 's, or equivalently, the corresponding active constraints in (|35]) 
with equality, via applying the property of polymatroid structure of CuAci'^*) given in (|34l) [22]. 

V. Simulation Results 

In this section, the performance of the proposed OMD is evaluated in comparison with the conventional 
SUD in a decentralized MU-MIMO system with K = 2 users, where the two users adopt an IWF- 
like algorithm to successively in turn optimize their transmit covariance matrices for individual rate 
maximization by deploying OMD or SUD at their receivers. For the purpose of exposition, all the 
channels involved in the system, including user's direct-link and cross-link channels, are assumed to 
have independent Rayleigh-fading distributions, i.e., each element of the channel matrix is independent 
and identically distributed as zero-mean CSCG random variable. Furthermore, each element of the two 
users' direct-link channels is assumed to have the variance pu and P22, for user 1 and 2, respectively; 
and each element of the two cross-link channels has the variance, pi2 for the channel from user 1 to 
user 2 and p2i for the channel from user 2 to user 1, respectively. In total, 5000 independent channel 
realizations are simulated over which each user's achievable average rate is computed. For each channel 
realization, the two users iteratively update their transmit covariance matrices until their rates both get 
converged. It is assumed that = = 2, k = 1,2. 

In Fig. m the achievable average sum-rate of the two users is shown for a symmetric system and 
channel setup, where Pi = P2 = 100, pu = P22 = 1, and pi2 = P21 = P- The user sum-rate is plotted 
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against p to investigate the effect of the interference between the two users on their achievable sum- 
rate. It is observed that the sum-rate with the proposed OMD improves over that with the conventional 
SUD for all the values of p, while the rate gains become more substantial in the case of large values 
of p, i.e., the "strong" interference case. With SUD, it is observed that the sum-rate first decreases with 
increasing of p (as a result of interference whitening), and then starts to increase with p (as a result of 
interference avoidance), and finally gets converged for large values of p (due to the fact that zero-forcing 
(ZF) -based receive beamforming to completely null the co-channel interference becomes optimal at the 
high signal-to-noise ratio (SNR) region). However, the sum-rate with the proposed OMD is observed to 
increase consistently with p, due to the fact that when the co-channel interference becomes stronger at 
the receiver, the OMD more easily decodes the interference. 

Next, we consider a special scenario of the general system model studied in this paper. In this case, 
a "cognitive radio (CR)" type of newly emerging wireless system is considered, where user 1 is the 
so-called primary (non-cognitive) user (PU) who is the legitimate user operating in the frequency band 
of interest, while user 2 is the secondary (cognitive) user (SU) that transmits simultaneously with the PU 
over the same spectrum under the constraint that its transmission will not cause the PU's transmission 
performance to an unacceptable level [23]. The PU is non-cognitive since it is oblivious to the existence 
of the SU and applies the conventional SUD at the receiver by treating the interference from the SU 
as additional noise. While for the SU, it is cognitive in the sense that it is aware of the PU and thus 
transmits with a much lower average power than that of the PU in order to protect the PU; thus, for this 
example it is assumed that Pi = lOP and P2 = P, where F is a given constant. In addition, since the SU 
is cognitive, it may choose to use the more advanced OMD at the receiver to cope with the interference 
from the PU. Two cases are thus studied for this example: Case (I) both user 1 and user 2 employ SUD; 
and Case (II) user 1 employs SUD while user 2 employs OMD. It is assumed that the SU's link distance 
is much shorter than that of the PU link, and furthermore the SU transmitter and receiver are both in the 
vicinity of the PU transmitter while they are both sufficiently far away from the PU receiver. Thus, for 
this example we assume that pu = l,p22 = 10,pi2 = 10, and P21 = 1- 

In Fig. [3l the achievable user individual rates are shown for different values of P in both Cases I and 
II. It is observed that the achievable rate of user 2 (the SU) improves significantly in Case II over Case 
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I, thanks to the use of OMD instead of SUD. This rate gain is substantial because the SU receiver is 
close to the PU transmitter and thus pi2 is large, i.e., the cross-link channel from PU to SU is a "strong" 
interference channel, for which the OMD is crucial for the SU to mitigate the PU's interference. However, 
it is also observed that the achievable rate of user 1 (the PU) drops slightly in Case II as compared with 
Case I. This is because that in Case II with OMD, the SU's transmitted signal has a more spatially 
spread-out spectrum than that in Case I with SUD, and so does the received SU's interference at the PU 
receiver. Nevertheless, due to the small value of p2i or the weak cross-link channel from SU to PU, the 
capacity loss of the PU is not significant, which justifies the operation principle of the SU, i.e., the PU 
transmission should be sufficiently protected. 

VI. Conclusion 

This paper studied a new decoding method, namely opportunistic multiuser detection (OMD), for 
the decentralized MU-MIMO system where each user iteratively optimizes transmit covariance matrix 
for individual rate maximization. In comparison with the conventional single-user detection (SUD), the 
proposed OMD still allows a fully decentralized processing of each user in the system, while it improves 
the user's interference mitigation capability at the receiver, and leads to more optimum spatial spectrum 
sharing among the users. Simulation results showed that substantial system throughput gains could be 
achieved by the proposed OMD over the conventional SUD, for certain application scenarios. 

Appendix I 
Proof of Lemma [3711 

We will prove Lemma [311 bv contradiction. First, suppose that /i* + /i2 < 1- Then, in the maximization 
problem of (fT6l) , from the expression of C{ri, Si, fii, [12) in (fTSl ). it follows that the optimal ri that 
maximizes the Lagrangian is = +00, which contradicts the fact that ri in (P2) is upper-bounded by 
finite rate values in the constraints (fTTl) and (fT2l) . Second, suppose that /ii + > 1- Similarly like the 
previous case, it can shown that rl = 0. However, this can not be true since we can easily find a feasible 
solution set for {ri,Si) in (P2) such that ri > 0. By combining the above two cases, it follows that 
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Appendix II 
Proof of Lemma [3T2] 



We rewrite R^^'^ in ^ and ^^"^ in ^ as 



= log 1/ + HnSi'^Hf, + H2iS2H^,\ - log 1/ + HuSi'^Hf, 



-f^ I 



?JD tx/f I 



= log 1/ + HnSrHf-i + H2iS2H^A - log 1/ + HnS^Hf^ 



I 



:.SD Tj-H 



(38) 
(39) 



Since 5^° and are optimal for the sum-capacity (in an equivalent two-user MIMO-MAC) and user' 1 
channel capacity (without the presence of user 2), respectively, we have 



log I J + I/fi + H2iS2Hg\ > log |/ + HuSfH^^ + H2iS2Hg\ 

log |/ + ifn^f J^nl < log |/ + HuSfH^il ■ 
Combining the above two inequalities with (|38l) and (l39l) . it thus follows that ^2"^ > R^\ 



(40) 
(41) 



Appendix III 
Proof of Proposition 14.11 



We first prove the former part of Proposition 14.11 i.e., the set is unique, by contradiction. Suppose 
that there exist two optimal decodable user sets for user 1 with the same size, denoted by Ai and 
Bi. Without loss of generality, we let Ai = {1^,0} and Bi = {£,C}^ where C, V and 8 are subsets 
consisting of completely different user indexes. Then, we can express Ai = {£,J^} and Bi = {V^J^}, 



where = Ai[JBi. Then, for users in the set Ai, their transmit rates must satisfy [20] 

-1 



^ Ti < loj 



keAi 



yj CV,}CCC. (42) 



Similarly, for users in the subset £ of Bi, their transmit rates must satisfy 



^ rj < log 



-1 



(43) 



Let J'' be an orthogonal set of J', where J^' [J = V. Similarly, I' is defined for X, where X' [Jl = £. 
(l42l) and (|43|) can thus be further shown as follows: 



^ Ti < log 



/+(/+ Yl H,,SkH^A Yl HaS.Hf, 

keilJT J ieJlJK 



(44) 
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^ Ti < log 



-1 



(45) 



From dSl) and (gS]), we obtain 



E 



n < lo£ 



+ log 



-log 



iexu^ 



log 



(46) 



Since 



log 



log 



HiiSiH^^ 



< lo£ 



lo£ 



(47) 



From (|46l) and (147]), it follows that 



iejnj/cijx 



r j < log 



log 



1+ J2 HaS.Hf, 



-log 



(48) 



Thus, the set |J /C |J X is a decodable user set for user 1 for any J' C V, JC C C, and ICE, and 
so is the set Qi = V[JC[J£. Since the size of Qi is larger than that of Ai or Bi, this contradicts the 
assumption that Ai and Bi are optimal decodable user sets for user 1. The proof of the former part of 
Proposition 14.11 thus follows. 

Next, we prove the latter part of Proposition 14. 1[ i.e., any decodable user set for user I, Ui, must be 
a subset of U*. The proof is also obtained via contradiction. Suppose that there is a set Ui that is not a 
subset of U^. Without loss of generality, we can express Ui = {V, C} and = {£, C}, where C, V and 
8 are orthogonal subsets. Based on the proof for the former part of Proposition 14.11 we know that the 
set P IJ C IJ £^ is also a decodable user set for user 1, and apparently, it has a larger size than til, which 
contradicts the fact that is the optimal decodable user set for user 1 . The proof of the latter part of 
Proposition 14.11 thus follows. 
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Initialize V = {2, . . 




While |V| > do 


(1) 


Initialize n = 1 




While n < 2^^! 


- 1 do 


If E»ev„ '^^ 


< C(V„) 


Set n ^ 


n + 1 


Else 




Set V ^ 


V- V„ 


Set 


VUVn 


Go to (1) 


End If 




End While 




Go to (2) 




End While 




Set Ul = V. 


(2) 



TABLE I 



The algorithm to find W. 



Appendix IV 
Algorithm to Find 

In this appendix, we present an algorithm to find the optimal decodable user set for user 1, Ul. First, 
some notations are given as follows for the convenience of presentation. Let Vn denote a subset of an 
arbitrary set V, n = 1, . . . , 21^' — 1. Note that here we have excluded the case that V„ = for the ease 
of presentation. The operation V — Vn then stands for removing the subset V„ from V. 

For a given user set, V C {2, . . . , K}, we know from Definition 14.11 that V is a decodable user set for 
user 1 if and only if for any subset of V, V„, it satisfies that 

5^ r, < log / + I / + ^ HkiSkH^, I Yl H^iS.Hf, ^ C{Vn). (49) 

However, if there exists a subset V„ such that X]jev„ ^« ^ C(Vn), it follows that V should not be a 
decodable user set for user 1. From the above property, we are able to design an iterative algorithm to 
find W*, which is explained as follows. Initially, we let V = {2, ... , K^. Thus, V = 0. Then, we will 
sequentially check for all the subsets of V whether Xljev„ ^« — C'(Vn),Vn. If this is the case, then we 
declare that = V. However, if we find any v! such that X^igv / '^^ ^ then we conclude that V 

should not be lAl and furthermore Ui CV — Vn'- In this case, we will set V ^ V — V„', V ^ V |J V„', 
and start a new sequence of tests for J2ieVn '"'^ — The above procedure iterates until we find 
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a set V such that YliieVn ^« — C'(V„), Vra or V = 0. In both cases, we set Ul = V. The above algorithm 
is summarized in Table HI 
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Fig. 1. The maximum achievable rate of user 1 with OMD, n, as a function of user 2's rate, r2, for some fixed 82. 
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2. The achievable sum-rate versus the average cross-link channel power gain, p, for a MU-MIMO system with K = 2, = 
= 1, 2, and Pi = P2 = 100. 




3. The achievable rate versus the average transmit power, P, in a MIMO CR system with Mk = Nk = 2, = 1, 2, Pi = lOP, and 
= P, for different decoding methods: Case (I) both PU and SU employ SUD; and Case (H) PU employs SUD and SU employs OMD. 



